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ON THE CRITICAL PARAMETER OF INTERLACEMENT 
PERCOLATION IN HIGH DIMENSION 

By Alain-Sol Sznitman 
ETH Zurich 

The vacant set of random interlacements on Z d , d > 3, has non- 
trivial percolative properties. It is known from Sznitman [Ann. Math. 
171 (2010) 2039-2087], Sidoravicius and Sznitman [Comm. Pure Appl. 
Math. 62 (2009) 831-858] that there is a nondegenerate critical value 
u* such that the vacant set at level u percolates when u < «» and 
does not percolate when u > it*. We derive here an asymptotic up- 
per bound on u,, as d goes to infinity, which complements the lower 
bound from Sznitman [Probab. Theory Related Fields, to appear]. Our 
main result shows that u* is equivalent to logd for large d and thus 
has the same principal asymptotic behavior as the critical parame- 
ter attached to random interlacements on 2d-regular trees, which has 
been explicitly computed in Teixeira [Electron. J. Probab. 14 (2009) 
1604-1627]. 

0. Introduction. Random interlacements have proven useful in under- 
standing how trajectories of random walks can create large separating inter- 
faces; see [4, 19, 20]. In the case of d > 3, it is known that the interlace- 
ment at level u > is a random subset of which is connected, ergodic 
under translations and infinite when u is positive; see [18]. The density of 
this set monotonically increases from to 1 as u goes from to oo. Its com- 
plement, the vacant set at level u, displays nontrivial percolative properties. 
There is a critical value U* in (0,oo) such that for u<u*, the vacant set at 
level u has an infinite connected component which is unique (see [16, 22]) 
and, for u > u*, only has finite connected components; see [18]. Little is 
known about and only recently was it shown that u* diverges when the 
dimension d tends to infinity; see [21]. The aim of the present article is to 
establish that u* is equivalent to logd as d tends to infinity. In particular, 
this result shows that has the same principal asymptotic behavior for 
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large d as the corresponding critical parameter (which has been explicitly 
computed in [23]) attached to the percolation of the vacant set of random 
interlacements on 2ci-regular trees. 

We now describe the model. Precise definitions and pointers to the lit- 
erature appear in Section 1. Random interlacements are made of a cloud 
of paths, which constitute a Poisson point process on the space of doubly 
infinite Z rf -valued trajectories modulo time shift, tending to infinity at pos- 
itive and negative infinite times. The nonnegative parameter u mentioned 
above plays the role (roughly speaking) of a multiplicative factor of the in- 
tensity measure of the Poisson point process. Actually, one simultaneously 
constructs, on a suitable probability space (fi,.4,P), the whole family I u , 
u > 0, of random interlacements at level u > [cf. (1.30)]. They are the 
traces on Z rf of the trajectories modulo time shift in the cloud having labels 
at most u. The complement V u of X u in Z d is the vacant set at level u. It 
satisfies the following identity: 

(0.1) P[V U 5 K] =exp{-ucap(iT)} for all finite K C Z d . 

In fact, this formula provides a characterization of the law on {0, l} zd of the 
indicator function of V u (cf. (2.16) of [18]). From Theorem 3.5 of [18] and 
Theorem 3.4 of [16], one knows that there is a critical value in (0, oo) 
such that: 

(i) for u > u*, P-a.s. all connected components of V u are finite; 
(0.2) (ii) for u <u*, P-a.s. there exists an infinite connected 



From Theorem 0.1 of [21], one has the following asymptotic lower bound on 
as d tends to infinity: 



The main aim of the present article is to show that the above lower bound 
does capture the correct asymptotic behavior of and that the following 
statement holds. 

Theorem 0.1. 



component in V". 



(0.3) 



lim inf u* / log d > 1 . 



d 



(0.4) 



lim u* / log d = 1 . 



As a byproduct, this result shows that it* has the same principal asymp- 
totic behavior as the critical value attached to random interlacements on 
2(f-regular trees when d goes to infinity; see Proposition 5.2 of [23]. We refer 
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the reader to Remark 4.1 for more on this matter. In addition, the proof of 
Theorem 0.1 also shows (cf. Remark 4.1) that 

(0.5) limu^/loga! = 1, 

d 

where u** G [u*,oo) is another critical value introduced in [19]. Informally, 
u** is the critical level above which there is a polynomial decay in L for the 
probability of existence of a vacant crossing between a box of side length L 
and the complement of a concentric box of double side length. It is an im- 
portant and presently unresolved question whether u* = u** actually holds. 
However, it is known that the connectivity function of the vacant set at level 
u, that is, the probability that and a (distant) x are linked by a path in 
V" (i.e., the probability of a vacant crossing at level u between and x) has 
a stretched exponential decay in x when u is bigger than it**; see Theorem 
0.1 of [17]. 

We will briefly comment on the proof of Theorem 0.1. In view of (0.3), 
we only need to show that 

(0.6) lim sup u* / log d < 1. 

d 

As for Bernoulli bond or site percolation, similarities between what happens 
on Z d and on 2<i-regular trees for large d lurk in the background of the 
proof. The statement corresponding to (0.6) for Bernoulli percolation is an 
asymptotic lower bound for the critical probability (a lower bound, not an 
upper bound, because the density of V u decreases with u), whereas the 
required lower bound in the Bernoulli percolation context follows from a 
short Peierls-type argument (cf. [3], page 640, [11], page 222, or [8], page 
25); the proof of (0.6) for random interlacements is quite involved. The 
long-range dependence present in the model is deeply felt. 

An important feature of working in high dimension is that the i 1 -, Eu- 
clidean and ^°°-distances all behave very differently on Z d ; see (1.1). At large 
enough scales (i.e., Euclidean distance at least d), the Green function of the 
simple random walk "feels the invariance principle" and is well controlled 
by expressions of the type (cyd/\ ■ \) d ~ 2 , where c does not depend on d 
and | • | stands for the Euclidean norm; see Lemma 1.1. However, at shorter 
range, the walk feels more of the tree-like nature of the space and the use 
of bounds involving the ^-distance becomes more pertinent [cf. (1.14) and 
Remark 1.2]. 

The above dichotomy permeates the proof of (0.6). We use a modifica- 
tion of the renormalization scheme ("for fixed d") employed in [17]. The 
renormalization scheme enables us to transform certain local controls on the 
probability of vacant crossings at level u$ = (1 + 5e) logci, e > 0, small, into 
controls on the probability of vacant crossings at arbitrary large scales at a 
bigger level Mqo < (1 + lOe) logtf. 
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The local estimates entering the initial step of the renormalization scheme 
are developed in Section 3. These involve controls on the existence of vacant 
crossings moving at -^-distance c(e)d from a box of side length Lq = d for 
the interlacement at level uq. The 2d-regular tree model lurks behind the 
control of these local crossings. The key estimates appear in Theorem 3.1 
and Corollary 3.4. These estimates result from an enhanced Peierls-type 
argument involving the consideration of what happens in -p £ 1 -balls, each 
having an £ 1 -radius c'ed and lying at mutual ^-distances of at least c"d. For 
this step, part of the difficulty stems from the fact that the local estimates 
need to be strong enough to overcome the combinatorial complexity involved 
in the selection of the dyadic trees entering the renormalization scheme. 

The renormalization scheme is developed in Section 2. It propagates along 
an increasing sequence of levels u n , with initial value uq = (1 + be) log d and 
limiting value Mqo < (1 + We) logd, uniform estimates on the probability of 
events involving the presence of certain vacant crossings at level u n . Roughly 
speaking, these events correspond to the presence in 2 n boxes of side length 
Lq(= d) of paths in V" n . The boxes can be thought of as the "bottom leaves" 
of a dyadic tree of depth n and are well "spread out" within a box of side 
length 3L n , where L n = IqLq and £o = d. The paths start in each of the 
2 n boxes of side length Lq and move at Euclidean (and hence i 1 -) distance 
of order c{e)d from the boxes. The estimates are conducted uniformly over 
the possible dyadic trees involved (cf. Propositions 2.1 and 2.3). The main 
induction step in the above procedure (cf. Proposition 2.1) relies on the 
sprinkling technique introduced in [18] to control the long-range interactions. 
The rough idea is to introduce more trajectories in the interlacement by 
letting the levels slightly increase along the convergent sequence u n . In this 
way, one dominates the long-range dependence induced by trajectories of 
the interlacement traveling between distant boxes. In the present context, 
the method uses, in an essential way, quantitative estimates on Harnack 
constants in large Euclidean balls when the dimension d goes to infinity. 
These estimates crucially enter the proof of Proposition 2.3. The bounds on 
the Harnack constants are derived in Proposition 1.3 with the help of the 
general Lemma A. 2 from the Appendix, which is an adaptation of Lemma 
10.2 of Grigoryan and Teles [7]. 

Let us now describe how this article is organized. 

In Section 1, we introduce notation and recall several useful facts concern- 
ing random walks and random interlacements. An important role is played 
by the Green function bounds (see Lemma 1.1) and by the bounds on Har- 
nack constants; see Proposition 1.3. 

In Section 2, we develop the renormalization scheme. It follows, with a 
number of changes, the general line of [17]. The key induction step appears 
in Proposition 2.1. The main consequences of the renormalization scheme 
for the proof of Theorem 0.1 are stated in Proposition 2.3. 
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In Section 3, we derive the crucial local control on the existence of vacant 
crossings at level uq traveling at ^-distance of order some suitable multiple 
of d. This local control is stated in Theorem 3.1. It enables one to produce 
the required estimate to initiate the renormalization scheme. This estimate 
can be found in Corollary 3.4. 

Section 4 provides the proof of (0.6). Combined with the lower bound 
(0.3) from [21], this yields Theorem 0.1. In Remark 4.1, we discuss some 
further questions concerning the asymptotic behavior of for large d. 

In the Appendix, we first derive, in Lemma A.l, an elementary inequality 
involved in proof of the Green function bounds from Lemma 1.1. We then 
present, in Lemma A. 2, a general result of independent interest providing 
controls on Harnack constants in terms of killed Green functions for general 
nearest-neighbor Markov chains on graphs. 

Finally, let us explain the convention we use concerning constants. Through- 
out the text, c or d denote positive constants with values which can change 
from place to place. These constants are independent of d. The numbered 
constants cq,c\, . . . are fixed as the values of their first appearances in the 
text. Dependence of constants on additional parameters appears in the no- 
tation, for instance, c(e) denotes a constant depending on e. 

1. Notation and random walk estimates. In this section, we introduce 
further notation and gather various useful estimates on simple random walk 
on Z d for large d. Controls on the Green function and on Harnack constants 
in Euclidean balls play an important role in the sequel. These can be found in 
Lemma 1.1 and Proposition 1.3. We also recall several useful facts concerning 
random interlacements. 

We let N = {0, 1,2,.. .} denote the set of natural numbers. Given a non- 
negative real number a, we let [a] denote the integer part of a. We denote 
by | • |i, | • | and | • |oo the £ 1 -, Euclidean and f°°-norms on M. d , respectively. 
We have the following inequalities: 

(1.1) |-|oo<|-|<|-|l, \-\<Vd\-\oo, \-\l<Vd\-\. 

Unless explicitly stated otherwise, we tacitly assume that d > 3. 

By finite path, we mean a sequence xq,...,xn in Z d , with N >1, which 
is such that — = 1 for < i < N. We sometimes write "path" in 
place of "finite path" when this causes no confusion. We denote by B(x,r) 
and S(x,r) the closed ball and the closed sphere, respectively, with radius 
r > and center x € In the case of the £ p -distance where p = 1 or oo, the 
corresponding objects are denoted by B p (x,r) and S p (x,r). For A,B C Z d , 
we write A + B for the set of x + y with x in A and y in B, and d(A, B) = 
inf{|x — y\;x G A, y G B} for the mutual Euclidean distance between A and 
B. We write d p (A,B), where p = 1 or oo, when the ^-distance is used 
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instead. The notation K CC Z d indicates that K is a finite subset of Z d . 
When £7 is a subset of Z d , we write |£7| for the cardinality of U, dU = 
{x G U c ;3y G £7, \x — y\i = 1} for the boundary of £7 and <9i nt £7 = {x G £7; 
By G £7 C , |x — 2/| i = 1} for the interior boundary of [7. We also write £7 in 
place of UUdU. 

We denote by W + the set of nearest-neighbor Z rf -valued trajectories de- 
fined for nonnegative times and tending to infinity. We write W+ and X n , 
n > 0, for the canonical cr-algebra and the canonical process on W+, re- 
spectively. We denote by 9 n , n > 0, the canonical shift on W + so that 
9n{w) = w(- + n) for u> G and n > 0. Since d > 3, the simple random walk 
on Z rf is transient and we write P x for the restriction to the set W+ of full 
measure of the canonical law of the walk starting at x G Z d . When p is a mea- 
sure on Z d , we denote by P p the measure ^xez d p{ x )Px and by E p the cor- 
responding expectation. Given U C Z d , we write iJjy = inf{n > 0;X n G £7}, 
= inf{n > 1; X n G £7} and Tjj = inf{n > 0; X n ^ £7} for the entrance time 
in £7, the hitting time of £7 and the exit time from £7, respectively. In the 
case of a singleton {x}, we simply write H x and H x for simplicity. 

We let g(-,-) stand for the Green function: 

(1.2) g{x,x) = ^P x [X n = x] forx,x'mZ d . 

n>0 

The Green function is symmetric in its two variables and, due to translation 
invariance, g(x,x') = g(x' — x) = g(x — x'), where 

(1.3) g(x) = g(x, 0) = 3(0, x) for x G Z d . 

The ^-distance is relevant for the description of the short-range behavior 
of g(-) in high dimension [cf. Remark 1.3(1) of [21] and Remark 1.2 below]; 
the Euclidean distance becomes relevant in the description of the "mid-to- 
long-range" behavior of <?(•). The following lemma will be repeatedly used 
in the sequel. We recall that the convention concerning constants is stated 
at the end of the Introduction. 

Lemma 1.1. 

(1.4) g(x) < (c Vd/\x\) d ~ 2 for\x\>d 
g(x) > ( Cl Vd/\x\) d ~ 2 

for \x\ 2 > d|x|oo > (and, in particular, when \x\ >d) 



(1.5) 

(1.6) P x [H B (fl,L) < 00] < ( o ) for L>d,x£Z d (with c>l) 



cL- d ~ 2 
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Proof. We begin with the proof of (1.4), (1.5). To this end, we de- 
note by pt(u,v), t > 0, u, v £ Z, the transition probability of the simple 
random walk in continuous time on Z with exponential jumps of param- 
eter 1. The transition probability of the simple random walk on Z rf with 
exponential jumps of parameter d can then be expressed as the product of 
one-dimensional transition probabilities. Relating the continuous- and the 
discrete-time random walks on Z d , we thus find that 

poo & 

(1.7) g{x) = d Y[pt(0,Xi)dt for x = ( Xl ,...,x d )£Z d . 
"' i=i 

From Theorem 3.5 of [15] and the fact that the function 

1 



F(rf) = - io g ( 7 + vV + 1) + -(vV + i-i), 7>0, 

7 

appearing in Theorem 3.5 of [15] has derivative —(1 + \/j 2 + l)" 1 , tends to 
in 7 = and thus satisfies the inequality log(l + 3r) < —F(j) < log(l + 7) 
for 7 > 0, we see that for suitable constants < n < 1 < k' , we have 

^(1 Vt V |u|)" 1/2 expj-Hlog^l + K -y^j I 

(1.8) <Pt{0,u) < ^(1 VtV|ii|)" 1/2 exp|-|u|log^l + K^-^| 

for t > 0, u G Z. 

We now prove (1.4) and thus assume that \x\ > d. By (1.7), (1.8), we bound 
g(x) from above as follows (we also use the inequality d <2 d and Lemma 
A.l from the Appendix): 

g(x) < ^^"(lVtr^expj-^l^llog^l + ^jj^ 
( < 1} c d J°°(l V t)- d / 2 exp|-|x| log(\ + I dt 



(1.9) 



< c d J^ X \l V t)~ d ' 2 expj - |x| log ^1 + J-^pj I dt 



+ c d fV d / 2 exp(-4^!>^ 



k\x\ 



■It 



where, in the last step, we have used the inequality log(l + 7) > \ for < 
7 < 1 . Performing the change of variable s = in the last integral, we see 
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that the last term of (1.9) is smaller than 

(1.10) c d \x\ 2 - d J ]X] 2 s d / 2 - 2 e- s ds<c d \x\ 2 - d v{^-l \ <(cVd/\x\) d - 2 , 

using the asymptotic behavior of the gamma function in the last step (cf. 
[14], page 88). 

As for the first integral in the last line of (1.9), we note that for 1 < s < 
k\x\, the function s — > — i logs — |sc| log(l + k^) has derivative 

d \x\ k\x\ s<k\x\ d U| \x\>d 

h— — -r-r > h— > 

2s s s + k\x\ 2s 2s 

and is hence nondecr easing. Thus, the first term in the last line of (1.9) is 
smaller than c d (K\x\)-^ 2 -^2-^ . 

Observe that for a > d, ^loga + a log 2 > (d - 2)log^ (indeed, this 

inequality holds for a = d and 4^ + log 2 > d ^ L for a > d) . It follows that 
the first term in the last line of (1.9) is at most {c^fd/\x\) d ~ 2 . Together with 

(1.10) , this completes the proof of (1.4). 

We now prove (1.5) and assume that x / 0. Since log(l + 7) < 7 for 7 > 0, 
and k' > 1, it follows from (1.7), (1.8) that 

g(x) > c d ^ t- d l 2 ^ V {-^\dt 

S = K'\x\ 2 /t Ax\ 2 /\X\ X 

(1.11) > c d \x\ 2 ~ d / s d / 2 ~ 2 e~ s ds 



, cVd\ d 2 i9 ,1 

> I -: — — I when \x\ > d\x\ 



\x\ 

and (1.5) follows. 

Finally, (1.6) is a routine consequence of the identity 

(1.12) g(x) = E x [g(XH B(0iL) ),H B{0>L) <oo} for L > and x G Z d , 
combined with (1.4), (1.5) and the fact that inf^( ,L) 9 — ^8B(o,L) 9- ^ 

Remark 1.2. 

(1) Although we will not need this fact in the sequel, let us mention that 
the following lower bound complementing (1.6) also holds: 

/ cL \ d ~ 2 

(1.13) Px[H B{0 ,L) < 00] > (j-rj A 1 for L > d, x G Z d (with c < 1). 

Indeed, one uses (1.12), together with (1.4), (1.5), and, when d+ 1 > L{> d), 
the inequality sup^ intB ( ,L) 9 — ^ sn VdB{o,L) 9i which follows from the fact 
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that g is harmonic outside the origin (the factor 2d can then be dominated 
by c d ~ 2 ). 

(2) Let us point out that when x = ([gP],0, ... ,0) with ^ < a < 1, the 
upper bound (1.4) does not hold when d > c(a). Indeed, it follows from 
(1.7), (1.8) that 

g(x) > dj\ t (0, [d a, ])pt(0 > 0) < *- 1 dt ^ c d d- a/2 exp{-d Q log(l + n'd 01 )}, 

which is much bigger than (coVd/\x\) d ~ 2 < c d_2 exp{ — (a — \){d — 2)logcfj- 
for d > c(a). 

(3) We recall from (1.11) of [21] that when d> 5, 

/ ,\d/2-2 

(1.14) g(x)<(j^) forxeZ d . 

The inequality is useful, for instance, when \x\ < d, but \x\i > C2d, a situation 
where (1.4) is of no help. We will use (1.14) in Section 3 when deriving local 
bounds on the connectivity function of random interlacements at a level uq 
close to logii; see the proof of Theorem 3.1. 

(4) The asymptotic behavior of g{x) for d fixed and large x is well known; 
see, for instance, [10], page 313, or [12], page 31: 

hm4^ = -r^-iV^ 2 . 



x-+oo \x\ d - 1 2 V 2 

The asymptotic behavior of g(-) at the origin, or close to the origin when 
d tends to infinity, is also well known; see, for instance, [13], page 246, or 
[21], Remark 1.3(1). On the other hand, the behavior of <?(•) at intermediate 
scales when d tends to infinity seems much less well explored. 

The bounds on the Green function of Lemma 1.1, together with Lemma 
A. 2 from the Appendix, enable us to derive quantitative controls on Har- 
nack constants in suitably large Euclidean balls. These bounds will be in- 
strumental for the renormalization scheme developed in the next section; see 
the proof of Lemma 2.2. First, we recall some terminology. When U C Z d , 
we say that a function u defined on U is harmonic in U if, for all x € U, 
u ( x ) = 2d S| e |=i u ( x + e )- We- can now state the following proposition. 

Proposition 1.3 (L > d). Setting c 3 = 4 + lOf 1 [where c > c\— see 
(l-4)> (1-5)], there exists c > 1 such that when u is a nonnegative function 
defined on B(0, C3L) and harmonic in B(0,C3L), we have 

(1.15) max u<c d min u. 

B(0,L) B(0,L) 
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PROOF. We define XJ X = B(0,L) QU 2 = B(0,4L) C U 3 = B(0, c 3 L). In 
view of Lemma A. 2 from the Appendix, any u as above satisfies the inequal- 
ity 



maxu < K minii, 

Ui Ui 



where 



(1.16) K= max max G U3 (x,z)/G U3 (y,z) 

and Gjj 3 (-,-) stands for the Green function of the walk killed outside U3 
[cf. (A. 8)]. Applying the strong Markov property at time Tu 3 and (1.2), we 
obtain the following identity: 

G Ua (y, z) = G(y, z) - E y [G(X Tu3 , z)] for y, z G Z d . 

Hence, when x,y & U\ and z G <9i n t^2, we see that 

( L4 > r- A 9 

(1.17) G Us (x,z)<G(x,z) < {c Vd/{2L)) d - 2 

and 

G U3 (y, z) > ( Cl Vd/(5L)) d - 2 - {c Vd/((c 3 - 4)L)}^ 2 



(1.18) 



d-2 / / \ d-2 / \ d-2 

Cl\ (Co 



L J \\5 J \C3-4 

f) (f) (-G 

We thus find that K < 2(^) d ~ 2 and the claim (1.15) follows. □ 

We now briefly review some notation and basic properties concerning the 
equilibrium measure and the capacity. Given K CC we write &k for the 
equilibrium measure of K and cap(i^) for its total mass, the capacity of K: 

e K (x)=P x [H K = oo]l K (x), x£Z d , 

(1.19) 

Cav (K) = Y / Px[H K = o }. 

The capacity is subadditive [a straightforward consequence of (1.19)]: 

(1.20) cap(KUK') < cap(K) + cap(K'), for K, K' CC Z d . 

One can also express the probability of entering K in the following well- 
known fashion: 

(1.21) P x [H K <^] = ^2g(x,y)e K (y) for x G 7* d . 

y&K 
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Further, we have the following bound on the capacity of Euclidean balls: 

(1.22) cap(5(0,L)) < f ^=J for L>d, 

which follows from (1.5), (1.6) and (1.21), by letting x tend to infinity. 

Remark 1.4. Although we will not need this estimate in the sequel, let 
us mention that in a way analogous with (1.4), (1.13) and (1.21), one finds 
that 



(1.23) cap(B(0,L)) > ( — 1 for L>d. 

We now turn to the description of random interlacements. We refer to 
Section 1 of [18] for details. We denote by W the space of doubly infinite 
nearest-neighbor Z rf -valued trajectories which tend to infinity at positive 
and negative infinite times and by W* the space of equivalence classes of 
trajectories in W modulo time shift. We let tt* stand for the canonical map 
from W into W* . We write W for the canonical u-algebra on W generated 
by the canonical coordinates X n , n G Z, and W* = {A C W*; (tt*)~ 1 (A) G 
W} for the largest <r-algebra on W* for which tt* : (W,W) — )■ (W*,W*) is 
measurable. The canonical probability space for random interlacements is 
now given as follows. 

We consider the space of point measures on W* x R + : 



n 

(1.24) 



< u = y^£( w y, Mi ),with (w*,Ui) G W* x R + ,i > and u > 0, 
^ i>0 



w(W* K x [0, u]) < oo for any K CC Z d and u > 



where, for K CC Z d , C W* stands for the set of trajectories modulo 
time shift that enter K, that is, = it*(Wk), where Wr is the subset of 
W of trajectories that enter K. 

We endow £1 with the <r-algebra A generated by the evaluation maps 
uj — > uj(D), where D runs over the u-algebra W* x Z3(R + ), and with the prob- 
ability P on (Q,A), which is the Poisson measure with intensity v(dui*)du 
giving finite mass to the sets x [0,u] for if CC u > 0, where v is 
the unique cr-finite measure on (VT*,W*) such that for any K CC Z d (see 
Theorem 1.1 of [18]), 

(1.25) I w ^u = tt* oQ K , 
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Qk here denoting the finite measure on W^, the subset of Wk of trajectories 
which are for the first time in K at time and such that for A, B E W+ [we 
recall that W+ is defined above (1.2)] and x E Z rf , 

Qk [(X_n.)n>0 E A, A" = X, (X n ) n > E 5] 

(1-26) 

= P a .[A|JT* = x]ex(:r)P a .[B]. 

For if EC Z d , u > 0, one defines on (Q,A) the following random variable 
valued in the set of finite point measures on (W+,W+): 

fiK, u {dw) = y~]8( w *)K,+ l{wZ E W* K ,Ui < u} 

i>0 

(1.27) 

for w = ^(5 (<illi) eO, 

i>0 

where, for E W^, (w*)^' + stands for the trajectory in W+ which follows 
w* step-by-step from the first time it enters K. 

When < v! < u, one defines fJ-K,u',u(dw) in an analogous way to (1.27), 
replacing the condition U{ < u with v! < Ui < u in the right-hand side of 

(1.27) . Then, for < u' < u, K EE Z d , one finds that 

f^K,u',u and fJ>K,u' are independent Poisson point pro- 

(1.28) cesses with respective intensity measures (u — u')P eK 
and u'P £K . 

In addition, one has the identity 

(1-29) HK,u = HK,u' + HK,u',u- 

Given lo E f2, the interlacement at level u > is the following subset of Z d : 
Z<» = (J ranged*) 

(1.30) 

if ^=EW)= U U 

i>0 KCCZ d w&Snppfj, KyU (uj) 

where, for w* E W*, range(w*) = w(N) for any w E W, with tt*(w) = w* , 
and Supp/ix,u(^) refers to the support of the point measure hk,u(^)- The 
vacant set at level u is the complement of T u {uj): 

(1.31) V' il (w) = Z d \Z M (u;) foruEO,u>0. 
One also has (cf. (1.54) of [18]) 

(1.32) X u (u)nK = [J w{N)nK for K C K' EE Z d ,u>0. 
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From (1.28), one readily finds that, as mentioned in (0.1), 

(1.33) P[V U 5 K] = exp{-ucap(iT)} for all K CC Z d , 

an identity that characterizes the law Q u on {0, l} zd of the indicator function 
of V u (ui); see also Remark 2.2(2) of [18]. This brings us to the conclusion of 
Section 1 and of this brief review of some useful facts that we will use in the 
following sections. 

2. From local to global: The renormalization scheme. In this section, we 
develop a renormalization scheme that follows, in its broad lines, the strategy 
of [17]. We introduce a geometrically increasing sequence of length scales 
L n , n > 0, and an increasing, but typically convergent, sequence of levels u n , 
n > 0. When the sequence u n is sufficiently increasing [cf. (2.19)], we are able 
to propagate from scale to scale bounds on the key quantities p n (u n ) that 
appear in (2.17). Roughly speaking, these controls provide uniform upper 
bounds on the probability that in a box at scale L n , 2 n "well spread" boxes 
at scale Lq all witness certain crossing events at Euclidean distance of order 
cLq in the vacant set at level u n . Interactions are handled by the sprinkling 
technique originally introduced in Section 3 of [18]. The renormalization 
scheme enables us to transform local estimates on the existence of vacant 
crossings at scale Lq in the vacant set at level uq into global estimates on 
crossings at arbitrary scales in the vacant set at level Uqo = limti n . The 
difficulty we encounter in the implementation of the scheme stems from the 
fact that we want both uo and Uoo to be "slightly above" the critical value 
(4.7) and (0.3). However, the local controls on vacant crossings at 
level uq, which we introduce into the renormalization scheme and develop 
in the next section, require Lq to be rather small, that is, of order d. We are 
then forced to keep a tight control on the estimates we derive when d goes to 
infinity. The Green function and entrance probability estimates from Lemma 
1.1, together with the bounds on Harnack constants in Euclidean balls from 
Proposition 1.3, play a pivotal role in this scheme. The fact that the £°°- and 
Euclidean distances behave very differently for large d [see (1.1)] also forces 
upon us some modifications of the geometric constructions in [17]; see, for 
instance, (2.1) and (2.26). The main results of this section are Proposition 
2.1, which contains the main induction step, and Proposition 2.3, which 
encapsulates the estimates we will use in Section 4. 

We consider the length scales 

(2.1) L >d, L = (Vd + R)L with R > 1 

as well as 

L n = PqLq 
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(2.2) for n > 1, where 4> > 1000— {\fd + R) is an integer 

c\ 

multiple of 100 (we recall that Co > c±; cf. Lemma 
1.1). 

We organize Z rf in a hierarchical way, with Lq being the finest scale and 
L\ < Li < • ■ • being coarser and coarser scales. Crossing events at the finest 
scale will involve the length scale Lq. We introduce the following set of labels 
of boxes at level n > 0: 

(2.3) I n = {n} x Z d . 

To each m = (n, i) E I n , n > 0, we attach the box 

(2.4) C m = (iL n + [0,L n ) d )nZ d . 
In addition, when n > 1, we define 

(2.5) C m = |J C ro /(DC ro ). 

m'6/«,d 00 (C m /,C m )<l 

On the other hand, when n = and m = (0, i) £ Io, we define instead 

(2.6) d m = B(iL ,L ) 1 ' 1 D 2 ' 1 (J S(x, J RL )(DC m ). 

xeC m 

The above definitions slightly differ from (2.3) in [17] due to the special 
treatment of the scale n = 0. It is relevant here to use Euclidean balls and, 
thanks to (1.6) of Lemma 1.1, to have a good control on the entrance prob- 
ability of a simple random walk in C m . The radius of these balls has to 
be chosen sufficiently large so that we can show that crossing events at the 
bottom scale, from C m to 9j n t C m , are unlikely (this will be done in the next 
section). 

We then write S m = d m tC m and S m = d m tC m for m G I n , n > 0. Given 
m G l n with 77, > 1, we consider 'Hi(m), Hiiyn) Q I n -i defined by 

Hi(m) = {m G I n -x; C m Q C m and CWnS m ^0}, 

(2.7) 

H 2 {m) = |m G J n _i; C w n L G Z d ; ^(z, C m ) = ^ | / 

We thus see that for n > 1, m £ I n , one has (see also Figure 1): 

m\£'Hi(m), m2 6%H 
(2-8) „ „ „ _ 

implies that C^j n Cm 2 = and CWi U Cm 2 C C m 

[in the case n = 1, we use the lower bound on £q in (2-2) as well as (1.1)]. 
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Fig. 1. An illustration of the boxes CW, and balls Gm i , i= 1,2, when m belongs to I\. 

Then, to each m G I n , n > 0, we associate a collection A m of "binary trees 
of depth n." More precisely, we define A m to be the collection of subsets T 
of Uo<fc<n Ik such that, writing T k = T n Ik, we have 

(2.9) T n = {m}, 

any m' G T" fe , 1 < k < n, has two "descendants," rfii(m') G Hi(m'), 



(2.10) 



i = 1,2, such that T k 1 = U m / 67 -fc{mi(m / ), m 2 (m')}. 



For each T G A m and m' G 7", one can then define the subtree of "descendants 
of m! in 7~" via 

(2.11) Tm! = {m" G T;C m » C C m ,}(GA m /). 

Given 1 < < n, m' G T fc , one thus has the following partition of T m i\ 

( 2 -12) 7^n' = U %n 1 (m') U 7^ 2 ( m /). 

In addition, we have the following rough bound on the collection A m of 
binary trees attached to m G I n : 

(2.13) |A m | < (c 4 4) 2(d - 1) (c 4 4) 4(d - 1) • • • (C44) 2 " (d - 1} = Mo) 2 ^^ 2 "- 1 ', 
where we have used the rough bound for m' G Ik, 1 <k <n, and, for i = 1,2, 



\Ui{m')\ < 2d[ c 



Lk 

L k -i 



d-l 



2d(ce ) d - 1 < {ciiof- 1 for some c 4 > 1. 



We then introduce, for u > 0, m£ I n , with n > 0, the event 



(2.14) 



^4 m — {Cm ^ ^ 'S'm}) 
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where the expression in the right-hand side of (2.14) denotes the collection 
of u) in £1 such that there is a path between C m and S m contained in V u . 
In an analogous fashion to Lemma 2.1 of [17], "cascades down to the 
bottom scale" because any path originating in C m and ending in S m must 
go through some Cfni 5 rn\ £ 1~L\{rfi} , reach Sjfi-^ and then go through some 
Cm 2 ,m2 £rl2(m), and reach SW 2 . Thus, similarly to Lemma 2.1 of [17], we 
find that defining for u > 0, n > 0, m 6 I n and T 6 A m 

(2.15) A\= p| A^, (recall that r° = rn/ ), 

m'eT 

one has the inclusion 

(2.16) AlQ (J Af. 

T eA ro 

We then introduce the key quantity 

(2.17) p n {u) = sup P[Aj-], u > 0,n > with m £ I n arbitrary, 

T eA m 

which is well defined due to translation invariance, and find that 

(2.18) F[A u m ]<\A m \ Pn {u) forn>0,n>0. 

The heart of the matter is now to find a recurrence relation bounding 
p n+ i(u n+ i) in terms of p n (u n ) for suitably increasing sequences u n (we are 
actually interested in increasing, but convergent, sequences). We recall that 
R appears in (2.1). 



Proposition 2.1. There exist positive constants 05,0$, c such that if 
£q > c(s/d + R), then, for any increasing sequences u n , n>0, in (0, 00) and 
nondecreasing sequences r n , n > 0, of positive integers such that 

( Ln (n+1)(d_2) \ r ™ +1 
(2.19) u n+1 >u n i 1 + J foralln>0, 



one has, for all n>0, 

p n+ l(u n+ i) <Pn(u n +l){pn(Un) + U n 

(2.20) 

X (Vl 



L 



Vd) 



(d-2) 



M ' r (n+l)(d-2) 
[note that p n (-) is nonincreasing so that p n (u n+ i) <p n {u n )]- 
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Proof. The proof of Proposition 2.1 is an adaptation of the proof of 
Proposition 2.2 of [17], which will be sketched below with some modifications 
which we will highlight. 

One considers some n > 0, m G I n +ii 7~ G A m and writes mi,m2 for the 
unique elements of Hi(m), H<i(jn) in T n {= T C\ I n ). One also writes u' 
and u, with < u' < u, in place of u n and u n +i. 

If T G A™ with m G I n , then one defines, for (i, a point process on W + 
defined on f2 (i.e., a measurable map from f2 into the space of point measures 
on W+): 



Ay(^) = [^] < w G f2; there is a path in 

C m / \ f [J io(N)^ joining C m / with S m '\. 



m'eT nl a 

(2.21) 



uiGSupp 

As in (2.19) of [17], using independence, we have the bound 

(2.22) P[A^]< Pn (u)F[A T2 (fi 2j2 )], 

where T2 stands for Tm 2 and we have decomposed the point process /iy. u 
[see (1.27)], where 

(2.23) V = Ci U C 2 

(2.24) with di= (J C m > C C mj for i = 1, 2 

(i.e., a union of 2" pairwise disjoint Euclidean balls of radius Lo)> into a sum 
of independent Poisson processes via 

(2.25) W,« = Mi,i +^1,2 + ^2,1 +A*2,2, 
where, for in {1,2}, we have set 

jUij = 1{X G Ci,HQ, < 00} fi v ,u and ^ = 1{A" G Ci^Hq^ = oo}fi v ,u- 

One introduces similar decompositions for fj,vu' m terms of analogously 
defined point processes /xj ■, 1 < i,j < 2, and for fiv,u',u in terms of //*^-, 
l<t,j<2. 

The heart of the matter is to bound P[A^= (^2,2)] = ^[^^(/^ 2 + ^2 2)]' 
which appears in the right-hand side of (2.22), in terms of p n (u') when u — u' 
is not too small. For this purpose, we employ the sprinkling technique of [18] 
and, loosely speaking, establish that ^ 2 dominates "up to small corrections" 
the contribution of ^x' 2 l + in PL4^ ] =P[A^ 2 (^ j2 + // 2 ,i + ^1,2)]- 
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With this in mind, we define a neighborhood U of Cfn 2 (and, in contrast 
to (2.20) of [17], we do not define U as the £°°-neighborhood of Cm 2 of size 
"t 1 ). Instead, if = {11,12) £ In [see (2.3)], we define U as the following 
Euclidean ball (which is much smaller than the corresponding ^°°-ball of 
same radius): 

(2.26) [/ = 5fi 2 L n ,^±ij DCW 2 using (2.1), (2.2), (2.5), (2.6). 
We then have the following important controls on Euclidean distances: 

dm c 2 ) > %tl - 3VdL n {2A) >' 2) %tl when n > 1 

(2.27) 10 ^° 
L n+1 ? (2.1) (2.2) L n+1 

> Ln > when n = 0, 

- 10 20 

and we have used in the first line the fact that C2 C Cm 2 when m > 1; see 
(2.8). Using similar considerations, we find that 

d(dU,C 1 )> ^±i-L n -L n -^±i-l>^±i whenn>l 

(2.28) 2 10 20 

> Ln — Ln 1 > when n = 0. 

~ 2 10 20 

Since V = C\ U C2, we have established that 

(2.29) dmV) >^. 

We then introduce the successive times of return to C2 and departure from 
U: 

(2.30) R\ = Hq , D\ = T\j o Rl + i? 1 and for > 1, by induction, 
R k+1 = R l oe Dk + D k , D k+1 = D 1 o6 Dk + D k , 

so that < Ri < Di < ■ ■ ■ < R k < D k < ■ ■ ■ < 00. 

Letting r > 1 play the role of r n in (2.19), (2.20), we further introduce 
the decompositions 

/4,1= ^2 P2,l+P2,l, tA,2= ^2 Pl,2 + Pl,2, 

l<i<r l<i<r 

(2.31) 

^2,2 = P2,2 + P2,2) 

\<Kr 
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where, for i ^ j in {1,2} and £>1, we have set 

Pi,j = l i R i <D e < Rt+i = oo}^, 

Pi,j = HRr+l <00}n' itj , 

pi 2 = 1{R(. <D e < R £+ i = Oo}/i 2>2 

and 

P2,2 = HRr+i < oo}p* 2 2 . 

The point processes p~i 2 and ~p 2 2 pl a Y the r °l e of correction terms, eventually 
responsible for the last term in the right-hand side of (2.20). The bounds we 
derive on the intensity measures £ 2 1 an d £1 2 °f P2 1 an d Pi 2 depart from 
(2.26), (2.27) in [17]. We write 

£2,1 (W+) = u'P ev [X e C 2 , < 00, Rr+i < 00] 

(1-19) 

(2.32) < u cap (C2) sup P x [Rr +1 < 00] 

Z6<?2 

strong Markov ^ / 

< ii / cap(C , 2)( sup Pa;^^ < OO 
Wee/ 2 

Combining (1.6) and (2.29), we find that 

(2.33) sup P x [H e2 < 00] < T (J*L\ ^ T (c^i^A ^ . 
Moreover, from (1.20), (1.22), we have 

(2.34) cap(a 2 )<2 n (^=J 
and hence 

/ f \ fd-2) / / f \ (d-2) \ r 

(2.35) < u' (-^ j ( 4 « (eg) j . 
In a similar fashion, we also obtain 

(2.36) ^(W + ) < u> (i* J ( 4 « (eg) ^(«+DC««> 

The next objective is to show that the trace on C2 of paths in the support of 
Si<£<r P21 an d Si<Kr /°i 2 * s stochastically dominated by the correspond- 
ing trace on C2 of paths in the support of /i 2 2 when ii — u' is not too small. 
An important step is the following lemma. 
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Lemma 2.2. For £ > c(\fd + R), all n>0, m £ 7 n+ i, T G A m , x € dU 
and y £ <9i nt C2, one /ias 

<i?i <oo,A^ 1 = y] 

(d-2) / \ (d-2)(n+l) 



(2.37) 



(2.38) 



P x [H^ < oo, Ri = oo 



f„\ (^-2) / _ \ (d-2)(n+l) 



Proof. The proof of (2.37) closely follows the proof of (2.30) in Lemma 
2.3 of [17]. The difference lies in the control of Harnack constants. Indeed, 
we first observe that the function h : z — > P Z [R\ < oo, X Rl = y] = P z [Hq^ < 

oo, Xh~ = y] is a nonnegative function, harmonic in Co. By (2.29), it is 

therefore harmonic on any B(zo, —^-) with zq G dU. One can then find 
c such that for any z, z 1 in dU, there exists a sequence Zi, < i < m, in 
dU with m < c, Zo = z , z m = z' and \zi + \ — Zi\ < in the notation of 

Proposition 1.3. Indeed, one simply projects z,z' onto the Euclidean sphere 
in R rf of radius L ^ 1 with center i2L n , the "center" of U [see (2.26)] and 
uses the great circle joining these two points to construct the sequence. 
Using (1.15) and a standard chaining argument, it follows that 

(2.39) sup P z [i?i < oo, X Rl = y]<c d inf P z [R 1<0 o, X Rl = y] . 

z&dU z ^ du 

The proof of (2.37) then proceeds as in Lemma 2.3 of [17] [and we use a 
similar bound to (2.33) above, where C\ replaces C2]. 

As for (2.38), we first note that for x G dU, due to (1.6) and (2.29), we 
have 

inf P x \Ri = 00, Hp, = ool 

x&dU Cl 

/ f \ ( d ~ 2 ) (2.2) / n 7 n \ ( d ~ 2 ) 

(2.40) >i_ 2 2"(c-M Vl-(f^ 

(2.1) 1 

> - v 

when £ >c'{Vd + R). 

On the other hand, a similar calculation leads to 



P x [H di < 00, Ri = 00] < 2 n ( c-"° 



->n+l 
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(2.41) 

, (Lo\ (d ~ 2) ( c\ {d - 2)(n+l) 

and (2.38) follows. □ 

The proof of Proposition 2.1 then proceeds like the proof of Proposition 
2.3 of [17] and yields that, under (2.19) (with v! in place of u n and u in 
place of u n+ i), 

(2.42) <p n {u')+l 2 ^W + )+l h2 (W + ) 

Inserting this inequality into (2.22), we thus infer (2.20) under the assump- 
tion on (2.19). □ 

We assume from now on that £q > c(y/d + R), with c > 2cs sufficiently 
large so that Proposition 2.1 holds. We then choose the sequences u n , n > 
and r n , n>0, as follows: 

(2.43) u n = u exp -M £ (r k + 1) if) 

^ ^ °' 0<fe<n ^ °' J 

(2.44) r n = r 2", 

where uq > and ro is a positive integer. The choice (2.43) ensures that 
(2.19) is fulfilled and the increasing sequence u n has the finite limit 

(2.45) u oo = „ e X p{(gH) < , ( 1 _ 2( J_ 1)( ,_ 2) + 1 _ (c ^ 1)( ,_ 2) )}. 

The next proposition reduces the task of bounding p n (u n ) to a set of con- 
ditions which enable us to initiate the induction procedure suggested by 
Proposition 2.1. We view a m as a function of uq, ro, £q, R [we introduced 
R in (2.1)]. 

Proposition 2.3. There exists a positive constant c such that when 
u > 0, r > 1, £ > c(Vd + R), L >d,L = (Vd + R)L , R>1 and K > 
log 2 satisfy 

/ f \ d-2 / q j \ r /2(d-2) 

(2 -> -(3) -MS) ■ 

(2.47) poM<e- Ko , 
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then 

(2.48) p n {u n ) < e -( A 'o-log2)2" for each n > Q 

Proof. The argument is similar to Proposition 2.5 of [17]. We assume, 
as mentioned before, that c > 2cs is large enough so that Proposition 2.1 
applies. Condition (2.46) implies that cqLq < £qLq(=Li). Thus, the last term 
in the right-hand side of (2.20) satisfies 

V '(.„( L y ' ( n+ i)( d _2) 



"•Us 4 ^' '» 



{d-2) / f- \ (d-2)r n / , \ nr. 



(2.49) <«oo -A c 6 



d-2 




(2.2), (2.46) / £ \ r»/2(d-2) 
< c 6 - 



As a result, (2.20) yields that for n > 0, 



£ r X (r /2)2"(d-2) 



(2.50) p n+ i{u n+ i) <Pn(u n )( p n (u n ) + ( C 6 -^- 

We then define by induction K n , n>0, via the following relation valid for 
n> 1: 

/ / / ?„ \ (ro/2)2»'(d-2)N 

(2.51) K n = K - £ 2-(™'+ 1 )log(l + e^ 2 " (c 6 /M ) 

so that X n < Kq and hence 

/ / / T n \ ( r o/2)2"'(d-2)x 

n'>0 V \ u u/ / 

(2.52) 

(2 > 6) if o - 2_(n ' +1) log 2 = K - log 2 > 0. 

n'>0 

As we now show by induction, we have p n (u n ) < e _i ^" 2 ". 

Indeed, this inequality holds for n = 0, due to (2.47), and if it holds for 
n > 0, then, due to (2.50), we find that 

/ / ? n \ (n)/2)2"-(d-2K 

P. +1 K +1 ) < (e— + (^) ) 



e -K„2»+i 1 + g ^ (c6 



^0 \ 

'toLo) 



(r /2)2 n (d-2) 



( 2 ^) e -K-„+i2»+i_ 
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This proves that p n (u n ) < e Kn2n for all n > and (2.48) follows. □ 

Remark 2.4. One of the main issues we now have to face is proving the 
local estimate Po(uq) < e~ K ° [see (2.47)] for large d, with uq of order close to 
log d (and a posteriori close to ) . We further need Kq sufficiently large so 
that 2 -( ^ 0-log2 ) 2n overcomes the combinatorial complexity arising from the 
choice of the binary trees in the upper bound (2.18), that is, overcomes the 

(2-13) 

factor | Aft | < (c4^o) 2 • Devising this local estimate will be the 

aim of the next section and will involve aspects of random interlacements 
at a shorter range, where features reminiscent of random interlacements on 
2d-regular trees (cf. Section 5 of [23]) will be evident. 

3. Local connectivity bounds. The aim of this section is to derive expo- 
nential bounds on the decay of the probability of existence of a path in the 
vacant set at level uq = (1 + 5e)logd, starting at the origin and traveling 
at ^-distance Md, where M is an arbitrary integer and d > c(e, M) (cf. 
Corollary 3.4). For this purpose, we develop an enhanced Peierls-type argu- 
ment. The main step appears in Theorem 3.1 below. In the present section, 
aspects of random interlacements on Z d for large d, reminiscent of random 
interlacements on 2d-regular trees (cf. [23]) will play a important role. We 
introduce the parameter 

(3.1) 0<e<i. 
We also introduce, in the notation of (1.14), 

(3.2) L = c 7 d where c 7 = [e 8 c 2 ] +2. 

The main result of this section is the following estimate on the connectivity 
function. 

Theorem 3.1 (d>c). For any positive integer M , we have 

v«o ( M(M - 1) e 2 

(3.3) P[0 «— > S\ (0, ML)} <exp<^ — i -L + 3Md Mdlogd 

[2 5 

where the notation is similar to (2.14) and 

(3.4) -u = (l + 5e)log(i. 

Proof. Observe that any self-avoiding path from to S\ (0, ML) succes- 
sively visits the ^-spheres <S*i(0, iL), i = 0, . . . ,M — 1. Thus, considering the 
first [jftd] steps of the path consecutive to the successive entrances in the var- 
ious spheres <Si(0, iL), we obtain M self-avoiding paths m, i = 0, . . . , M — 1, 
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where 7T, starts in Si(0,iL) and has [j^d] steps for each i. Denoting by Zi, 
i = 0, . . . , M — 1, the respective starting points of these paths, we find that 

(3.5) P[0^5i(0,ML)] < ^P[V U() D range vr* for i = 0, . . . , M - 1], 

where the above sum runs over Zi £ S(0,iL) and self-avoiding paths 7Tj with 
[j(jd] steps and starting points Zj, i = 0, . . . , M — 1. The next lemma provides 
a very rough bound on the cardinality of ^-spheres and ^ 1 -balls. Crucially, 
it shows that ^-spheres and balls of radius cd are "rather small," that is, 
their cardinality grows at most geometrically in d. 

Lemma 3.2 (£ e N). 

(i) \S 1 (0,£)\<2 d e e+d , 

(3.6) 

(ii) 151(0,01 <2 d e e+l+d . 

Proof. We express the generating function of |<Si(0, fc)|, k > 0, as fol- 
lows. Given Itl < 1, we have 



fc>0 fc>0 mi,...,m d >0 

m\-\ \-m^=k 



,d};m;^0}| 



(3.7) 



t mi ^ hm d 2l{«6{l,...,^};mi7^0}| 

mi,...,m,j>0 



m>l ' 



l + 2° 

< 



1-t J ~ (1 -t) d ' 

As a result, we see that for < t < 1, I > 0, 

|Si(0,£)| <2 d (i-t)- d r^. 
Choosing t = ^/(d + £), we find that 

(3.8) |Si (0,£)| <2 d (l + -\ + <2 d e e+d , 



dj V «. 

where we have used the inequality 1 + u < e u in the last step. This proves 
(3.6) (i). As for the inequality (3.6) (ii), by (3.6) (i), we can write 

I £+1 _ i 

(3.9) |Si (0,£) I < 2 d e dS Te k = 2 d e d - < 2 d e e+1+d 

e — 1 

fc=o 
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and our claims follows. □ 

We now come back to (3.5). By a very rough counting argument for the 
number of possible choices of 7Tj, we have a Peierls-type bound: 

ASi(0,ML)] 

< ( J] \s 1 (o,kL)\(2dy/ 10 A 

V fc=0 / 

x sup F[V U0 2 range m, i = 0, . . . ,M - 1} 



(3.6)(i) (M-l \ 
(3.10) < [ 2 d e kL+d 



{2df' WMd 

\ k=0 ) 

x supP[V"° D ranges, i = 0,...,M-l] 



< e Af(Af-l)/2L+2Afd/2^e/10Afd 

x supP[V"° D ranges, i = 0,...,M- 1], 

where the supremum runs over a similar collection as the sum in (3.5). 

The next objective is to bound the probability in the last line of (3.10). 
For this purpose, for each x in the set 

M-i 



B = U 



i=0 
(3.11) 

(pairwise disjoint ^ 1 -balls appear in this union), 

we write z x for the unique Zi such that x G -B^i, to^). We then define, for 
any x in B, the subset W* of PF* — see above (1.24) — (not to be confused 
with W? y): 

W* = the image under ir* of 

w 6 W: the minimum of di(z x ,w(n)), n 6 Z, 

is reached for the first time at u>(n) = x and w; 
does not enter any B\ ( Zi, -^d J with zi ^ z x 



(3.12) 



Note that, clearly, W* C and that 

(3.13) W*, x £ B, are pairwise disjoint measurable subsets of W* . 
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It then follows that for Zi,7Ti, < i < M — 1, as in (3.10), we have 
P[V M0 D range 7r is i = 0, . . . , M - 1] 

/Af-l \ 



< 



(3.14) 



[J U w;x[o,«o] =o 

i=0 xerangeTi / 



{Af-l 
-^E E kw:) 

f ed 
< expj -ti M- x inf u(W*' 

We will now seek a lower bound on v{W*) for x £ B. 

Choosing = {x} in (1.25), (1.26), by (1.19), we see that for any x in B, 

u{W;) = P x [\X n - z x \x > \x - for n > 0, 

and #U^fii(^/i0d) =O0 ] 
x P x [|A" n - 2^ 1 1 > |x - z x \i, for n > 0, 

> ( P x [|^n -z x \i>\x- z x \i, for n > 0] 



(3.15) 



E i *[ fl Bi(a tl e/10(0 <°°]) 



In view of (1.14) and the choice of L in (3.2), we see that when d > 8, we 
have 



E 



< 



(3.16) 



(1.14), (3.6)(ii) 

< 



E SU P 9 {y - x) B x ( 0, ^-d 

tf z y£B( Zl ,e/10d) V iU 



c 2 d \ 



2 Y<\jL-e/hd) 



d/2-2 



2d e e/10d+l+d 



J>1 

( < 2) 2 e - 8 ( d / 2 - 2 )+ M+1 Ei" (d/2 " 2) ? ce_d - 

i>i 

The next lemma yields a lower bound on the first term in the last line of 
(3.15). 
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Lemma 3.3 (d> c). When \y\i < | 3 one has 
(3.17) P y [\X n \ 1 >\y\ 1 foralln>0]>l 



2d-(|s/|i VI)" 



Proof. We first note that for z = (^i,...,«d) in Z d , P z -&.s., \\Xi\i 
\z\i\ = 1, and 



/ d 
|X 1 | 1 = |^| 1 + 1] = — 2d- 



(3.18) 



2d V ^ 

\ k=l 

where p m = f ( - + - [ 




We then introduce the canonical Markov chain N n on N that jumps tom+1 
with probability p m and to m-1 with probability q m = l— p m when located 
at m. We denote by Q m the canonical law of this Markov chain starting in 
to. In view of (3.18), a coupling argument shows that we can construct X n 
and N n on the same probability space so that a.s. \X n \x > N n for all n > 
and Xq = y £ lj d , No = \y\i. Consequently, we see that when y 7^ 0, we have 
the bound (with m = \y\\ < f ) 

PyWs^cP) < #Bi(0,|j/|i)] 
(3.19) >Q lyh [H d2 <H lyll ] 

= Pm(l + Pm+1 + Pm+lPm+2 H 1" Pm+1 ' " ' P<Z 2 -l) _ , 

where = ^ for £ > and we have used [5], (5), page 73. 

Note that the expression in the right-hand side of (3.19) is a decreasing 
function of each p#, m + 1 < £ < d 2 . If we further observe that pi < (i — \ 



x 



3X5 + 5 x 3) = f for m + 1 < ^ < f d and p £ < 1 for f d < £ < d 2 - 1, then 
we see that the above expression is bigger than 

7 v fc>0 v 7 v 7 



^ m\/ 5 m 5/3y /4 
2dyV + 2 2d-m + 3V5y 

-1 



d>c 



l<m<d/2 



fn \ 1 _ „ m 



> l-^ 1 + 3 



2d 7 V 2d - m 



m\ ( m \ ( d N 

> 1 1 - 3— >0 since m < - 

2d A 2d-m \~ ~ 2, 
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By the strong Markov property at time H Sl (o t d 2 ), we thus find that for d> c, 
1 < Ml < o> we have 

-Pj/[|^n|i > \y\i for all n > 0] 

(3-20) > (l - (l - ^jL-) (l - |z sup 2 P z [^ l(0jd/2) < oo] 

^( 1 -^)( 1 -^)- w ^^<°°l 

(1 > 6) i _ \vh _ 3 i^ii + M _ ^V d_2) 

2o! 2d-|y|i 2d(2d-|i/|i) Vv^7 



2d-|y| 



This completes the proof of (3.17) when y ^ 0. The extension to the case 
y = is immediate. □ 

We use the above lemma to bound the first term in the last line of (3.15) 
from below. In view of (3.16) and (3.17), we thus find that for d>c and any 
x£B [see (3.11)], 

(3.21) v(w;)>(i-5 'v^y 1 v>i-io-^l>i- £ . 

v ; k xJ-y 2d-(\x-z x \ 1 yl)J ~ 2-e/lQ- 
Coming back to (3.14), we thus find that 

(3.22) F[V U ° D ranges, i = 0,...,M-l] < exp|-^Me(l - e)d 

Inserting this bound into the last line of (3.10), we obtain 
P[0<^Si(0,ML)] 

< exp | M(M~l) L + _ ^ £(1 _ £)Md y 2d) e/10M d 

(3-4) (M(M-l) e 

< exp< — y — J -L + 3Md+ — Mdlogd 

Zi J.U 

__L( e + 4e 2 -5e 3 )Af(ilogd|. 

Since 5e 3 < 2e 2 , due to (3.1), the claim (3.3) follows. 

We will use the following corollary in the proof of Theorem 0.1 in the next 
section. 
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Corollary 3.4 [with (3.1), (3.2)]. If M>1, then for d > c{M , e) , 
(3.23) P[ott<Si(0,ML)] <exp|-^dMlogd|. 

Proof. This is an immediate consequence of (3.3). □ 

Remark 3.5. One should note that the bound of Theorem 3.1 deterio- 
rates when M becomes large. One can view Theorem 3.1 as a Peierls-type 
bound (slightly enhanced due to the role of M in the proof). In the next sec- 
tion, we will choose M as a large constant depending on e and use Corollary 
3.4 to produce the local estimate which will enable us to initiate the renor- 
malization scheme of Section 2. In this way, the local estimate on crossings 
in V M ° at ^-distance of order c(e)d will be transformed into an estimate on 
crossings at all scales in V w °° , where u^q <(l + 10e)logd. 

4. Denouement. In this section, we prove Theorem 0.1. We combine the 
local bound on the connectivity function at level uq of the last section (cf. 
Corollary 3.4) with the renormalization scheme of Section 2 (cf. Proposition 
2.3) in order to produce a bound on vacant crossings at a level Uoc £ [(1 + 
5e)logd, (1 + 10e)logd], valid at arbitrarily large scales. 

Proof of Theorem 0.1. We choose e and no as in (3.1), (3.4), respec- 
tively. For the renormalization scheme of Section 2, we choose [the constant 
C7 appears in (3.2)] 

(4.1) L = d, L = (Vd + R)L with R = 300c 7 e~ 2 
and 

(4.2) £ = d. 

In the notation of Proposition 2.3 and (2.13), we choose 

(4.3) r = 24 
and 

(4.4) K = log(4(c 4 4) 2(d - 1) ) ( = } log(4( C4 d) 2 ^ 1 )). 
In the application of Corollary 3.4, we choose 

(4.5) M = [100e~ 2 ] + 1 
so that in the notation of (3.2), (4.1), 

(4.6) ML + 1<RL . 
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We will now check that the assumptions of Proposition 2.3 hold for d > c(e). 
By (2.45), we see that for d > c(e), 

(4.7) u = (1 + 5e) log d < Uoo < (1 + lOe) log d 
and also that 

(4.8) L < 2d 3/2 . 
As a result, we find that 

(4.9) Uoo {^T ^ ^( 1 + 10e )( 1 °g d )( M ) (d " 2) 
and that 

(4.10) e^° = 4(c 4 d) 2(d - 1) , 
whereas, on the other hand, 

(4.11) >(cd) 6 ( rf - 2 ). 

Since 2{d— 1) < 6(d — 2), we see that for d > c(e), the expression in the left- 
hand side of (4.11) dominates the corresponding expressions in (4.9) and 
(4.10), that is, (2.46) holds. 

There remains to check (2.47). For this purpose, we apply Corollary 3.4 
and find that for d > c(e), since L > s/dL Q + ML + 1 [cf. (4.1), (4.6)], we 
have 

Po (u ) = P[[0,L - l] d d int B(0, Lq)] 

u (3 23) 

(4.12) <L^P[0^45i(0,ML)l < exp{d log d - lOci log cZ} 

(4.5) 

= (T M . 

We thus find that for d > c(e), po( u o) < e - ^ , that is, (2.47) holds as well. 
It now follows from Proposition 2.3 that for d > c(e), 

(4.13) Pn(uoc) < e-( K °- lo ^ 2n for all n > 0. 
Taking (2.13), (2.18) into account yields that for all n > 1, 

P[[0, L n - l] d f—> & mt [-L n , 2L n - l] d ] 

(4.14) 

< (C4£o) 2{d ~ 1){2n ~ 1) e"^°- los2 ) 2n (4 < 0) 2~ 2 ". 

In particular, the above inequality implies that P[0 < — > oo] = and hence 
u* < Uoo < (1 + 10e) logd for d > c(e). The claim (0.6) readily follows. Com- 
bining this upper bound with the lower bound (0.3), we have thus proven 
Theorem 0.1. □ 
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Remark 4.1. 

(1) The inequality (4.14), together with the fact that L n = Lq£q for n > 0, 
is more than enough to show that for e as in (3.1) and d > c(e), 

V (l+10e)logd 

lim mP (o, L) < — > Soo (0, 2L)] = 0, 

for some, and, in fact, all, 7 > 0. From the definition of the critical parameter 
■u** in [19], 

u** = inf{-u > 0; a(u) > 0} 

(4.15) 

wherea(u)=supfa>0; lim L a P[B 0C (0, L) Soo (0, 2L)} = oj 

(the supremum is, by convention, equal to zero when the set is empty), we 
thus find that for d > c(e) , 

(4.16) u**<(l + 10e)logd 

Since < u** , it follows that we have also proven that 

(4.17) limu**/logd= 1. 

d 

It is presently an open question whether u* = u**; however, we know that 
< < < 00 for all d > 3 (cf. [21]) and that for u > u^, the connectivity 
function has a stretched exponential decay (cf. [17]). 

(2) One may wonder whether the following reinforcement of (0.4) actually 
holds: 

P[0GV M *]=e-"* /3(0) ~(2d)- 1 asd^oo. 

This would indicate a similar high-dimensional behavior as for Bernoulli 
percolation; see [1, 2, 6, 9, 11]. In the case of interlacement percolation on a 
2ci-regular tree, such an asymptotic behavior is known to hold (cf. [22]). 



APPENDIX 

In this appendix, we prove an elementary inequality which is involved 
in the proof of the Green function estimate (1.14); see Lemma A.l below. 
We then prove, in Lemma A. 2, a bound on Harnack constants in terms of 
killed Green functions for nearest-neighbor Markov chains on graphs. The 
result is stated in a rather general formulation due to the fact that it is of 
independent interest. It is an adaptation of Lemma 10.2 of [7]. We recall 
that Lemma A. 2 enters the proof of Proposition 1.3. 
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Lemma A.l. 



Fora,b>0 Va 2 + b 2 log(l + \J a 2 + b 2 ) 

(A.l) 

<alog(l + a) +61og(l + 6). 

Proof. We introduce ip(u) = ulog(l + u), u > 0, as well as ipb( a ) = 
\J a 2 + b 2 and Xb( a ) = VK a ) + V'(^) — VKVbC )) for o, 6 > 0. We want to show 
that 

(A.2) Xb(a)>0 fora,6>0. 

We note that Xfe(0) = and that 

X ' b (a) = log(l + a) + 1 - — j— - (log(l + ^(a)) + 1 - 1 ) 

1 + a V 1 + <M a )/ 

The claim (A.2) will follow once we show that 

(A.3) X b (a)>0 fora,6>0. 

To this end, we note that for a > 0, x'o( a ) = an d that 

d . , . / 1 6 1 ft \ a 

-Xft(a) = - — TT^T + 



56 \l + tp b (a) (pb(a) (1 + ^fe(a)) 2 ^b(a) J p b (a) 

f, /-, , „ , 1 \ ab 

+ log(l + y> & (a)) + 1 



(A.4) 



1 + (^b(a)y (^(a) 3 
aft 



(1 + (pb(a))p b (a) 3 
x(log(l + V6 (o))(l + ^(a)) 



1 + 99 fe (a) 

We introduce the function p(u) = log(l + u)(l + u) — u>0. Observe that 
p(0) = and p'(tt) = log(l + u) + 1 - {l ^ u)2 > so that p(u) > for u > 0. 

Coming back to the last line of (A.4), we find that for a > 0, JjXf>( a ) — f° r 
ft > 0. This shows (A.3) and the claim (A.l) then follows. □ 

We now turn to the second result of this appendix. We consider a con- 
nected graph T with an at most countable vertex set E and edge set £ (a 
subset of the collection of unordered pairs of E). Given U C E, we define dU, 
di n tU and U similarly to what is described at the beginning of the Section 
1 (with obvious modifications). We consider an irreducible Markov chain on 
E, nearest-neighbor in the broad sense (i.e., at each step, the Markov chain 
moves to a vertex which is at graph-distance at most 1 from its current loca- 
tion). We write X n , n > 0, for the canonical process, P x for the canonical law 



CRITICAL PARAMETER OF INTERLACEMENT PERCOLATION 33 



starting from i£B and otherwise use similar notation as described at the 
beginning of Section 1. We denote by p{x,y), x,y 6 E, the transition prob- 
ability. We assume that the Markov chain satisfies the following ellipticity 
condition: 

(A. 5) p(x,y) > when x,y are neighbors (i.e., {x,y} 6 £). 
For / a bounded function on E, we define 

(A.6) Lf(x) = E x [f(X 1 )]-f(x) = Y,P(x,y)(f(y)-f(x)) iovxeE, 

where y ~ x means that y = x or y is a neighbor of x. Given U C E, a 
bounded function on U is said to be harmonic in U when (with a slight 
abuse of notation) 

(A.7) Lf(x) = for x € U. 

When U is a finite strict subset of E, the Green function killed outside U is 
defined as follows (the notation is similar to that in Section 1): 



(A.8) Gu(x,y) = E x 



J2Hx k =y,Tu>k} 



k>0 



x,yeE. 



It follows from the ellipticity assumption (A. 5) that when U is connected, 
Gjj(x,y) > for all x,y £U. The next lemma is an adaptation of Lemma 
10.2 of [7]. 

Lemma A. 2. Assume that ^ U\ C U2 Q Us are finite strict subsets of 
E, with Us connected, and that u is a bounded nonnegative function on Us 
which is harmonic in Us- We then have 

(A.9) maxu < A'min-u, 

f/i ~ I7i 

where 

(A. 10) K= max max Gu 3 (x,z)/Gu 3 {y,z). 

x,yeUi ze9i n tC/2 

PROOF. We define, for x <E E, 
(A.11) v(x) = E x [u{X Hu2 ),H U2 < T U3 ]. 

We first note that 

(A. 12) u(x)>v(x) foixGUs and u(x) = v(x) forx£U2- 
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Indeed, in view of (A. 11), u and v agree on U2 and, thanks to our assump- 
tions, u(X nA T U:j ), n > 0, is a bounded martingale under P x , x G U3, so that 
by the stopping theorem, we find that 

u(x) = E x [u(X Hu2ATu3 )] = v(x) + E x [u(X Tus ),T Us < H U2 ] 
> v{x) for x G U3. 

The claim (A. 12) then follows. 

Applying the simple Markov property at time 1 in (A.ll), when x G 
U3 \ U2, we see that 

(A. 13) v is harmonic in U-s\U2- 

In addition, we have, for x G U2, 

(A. 12) 

v(x) = u(x) = 22p(x,y)u(y) > 22p(x,y)v(y) 

y^x yr^x 

and the last inequality is an equality when x G U2 \ d- m t U2 ■ We have thus 
shown that 

(A.14) Lv = l dintU2 Lv<0 onU 3 . 

Applying the stopping theorem, we see that, under any P x , 

v(X nA T U3 ) — Lv{Xk), n > 0, is a martingale. 

0<fc<nAT C / 3 

Taking expectations and letting n tend to infinity, we obtain the identity 



v(x) = E x [v(X Tu3 )]-E : 



0<k<T U3 



(A.15) = -J2 G u 3 (x,z)Lv(z) 

(A i 4) J2 G U3 (x,z)(-Lv)(z), xeE. 

z&di nt U2 

Since v and u agree on C/2 2 (A.9) is a direct consequence of the above 
representation formula for v. □ 
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